Multi-scale statistics of turbulence motorized by active matter

A number of micro-scale biological flows are characterized by spatio-temporal chaos. These include dense suspensions of swimming bacteria, microtubule bundles driven by motor proteins, and dividing and migrating confluent layers of cells. A characteristic common to all of these systems is that they are laden with active matter, which transforms free energy in the fluid into kinetic energy. Because of collective effects, the active matter induces multi-scale flow motions that bear strong visual resemblance to turbulence. In this study, multi-scale statistical tools are employed to analyze direct numerical simulations (DNS) of periodic two- (2D) and three-dimensional (3D) active flows and compare them to classic turbulent flows. Statistical descriptions of the flows and their variations with activity levels are provided in physical and spectral spaces. A scale-dependent intermittency analysis is performed using wavelets. The results demonstrate fundamental differences between active and high-Reynolds number turbulence; for instance, the intermittency is smaller and less energetic in active flows, and the work of the active stress is spectrally exerted near the integral scales and dissipated mostly locally by viscosity, with convection playing a minor role in momentum transport across scales.Comment: Accepted in Journal of Fluid Mechanics (2017

Theory of the propagation dynamics of spiral edges of diffusion flames in von Kármán swirling flows

This analysis addresses the propagation of spiral edge flames found in von Kármán swirling flows induced in rotating porous-disk burners. In this configuration, a porous disk is spun at a constant angular velocity in an otherwise quiescent oxidizing atmosphere. Gaseous methane is injected through the disk pores and burns in a flat diffusion flame adjacent to the disk. Among other flame patterns experimentally found, a stable, rotating spiral flame is observed for sufficiently large rotation velocities and small fuel flow rates as a result of partial extinction of the underlying diffusion flame. The tip of the spiral can undergo a steady rotation for sufficiently large rotational velocities or small fuel flow rates, whereas a meandering tip in an epicycloidal trajectory is observed for smaller rotational velocities and larger fuel flow rates. A formulation of this problem is presented in the equidiffusional and thermodiffusive limits within the framework of one-step chemistry with large activation energies. Edge-flame propagation regimes are obtained by scaling analyses of the conservation equations and exemplified by numerical simulations of straight two-dimensional edge flames near a cold porous wall, for which lateral heat losses to the disk and large strains induce extinction of the trailing diffusion flame but are relatively unimportant in the front region, consistent with the existence of the cooling tail found in the experiments. The propagation dynamics of a steadily rotating spiral edge is studied in the large-core limit, for which the characteristic Markstein length is much smaller than the distance from the center at which the spiral tip is anchored. An asymptotic description of the edge tangential structure is obtained, spiral edge shapes are calculated, and an expression is found that relates the spiral rotational velocity to the rest of the parameters. A quasiestatic stability analysis of the edge shows that the edge curvature at extinction in the tip region is responsible for the stable tip anchoring at the core radius. Finally, experimental results are analyzed, and theoretical predictions are tested

Preface from the Guest Editors

This article is an editorial: Čantrak Đ., Urzay J. (2017/2018): Preface from the Guest Editors, Thermal Science, Vol. 21, Suppl. 3, pp. SVII, ISSN 2334-7163 (online edition), ISSN 0354-9836 (printed edition)

Diffusion-flame ignition by shock-wave impingement on a supersonic mixing layer

Ignition in a supersonic mixing layer interacting with an oblique shock wave is investigated analytically and numerically under conditions such that the post-shock flow remains supersonic. The study requires consideration of the structure of the post-shock ignition kernel that is found to exist around the point of maximum temperature, which may be located either near the edge of the mixing layer or in its interior, depending on the profiles of the fuel concentration, temperature and Mach number across the mixing layer. The ignition kernel displays a balance between the rates of chemical reaction and of post-shock flow expansion, including the acoustic interactions of the chemical heat release with the shock wave, leading to increased front curvature. The analysis, which adopts a one-step chemistry model with large activation energy, indicates that ignition develops as a fold bifurcation, the turning point in the diagram of the peak perturbation induced by the chemical reaction as a function of the Damköhler number providing the critical conditions for ignition. While an explicit formula for the critical Damköhler number for ignition is derived when ignition occurs in the interior of the mixing layer, under which condition the ignition kernel is narrow in the streamwise direction, numerical integration is required for determining ignition when it occurs at the edge, under which condition the kernel is no longer slender. Subsequent to ignition, for the Arrhenius chemistry addressed, the lead shock will rapidly be transformed into a thin detonation on the fuel side of the ignition kernel, and, under suitable conditions, a deflagration may extend far downstream, along with the diffusion flame that must separate the rich and lean reaction products. The results can be helpful in describing supersonic combustion for high-speed propulsion

Thermochemical effects on hypersonic shock waves interacting with weak turbulence

The interaction between a weakly turbulent free stream and a hypersonic shock wave is investigated theoretically by using linear interaction analysis (LIA). The formulation is developed in the limit in which the thickness of the thermochemical nonequilibrium region downstream of the shock, where relaxation toward vibrational and chemical equilibrium occurs, is assumed to be much smaller than the characteristic size of the shock wrinkles caused by turbulence. Modified Rankine-Hugoniot jump conditions that account for dissociation and vibrational excitation are derived and employed in a Fourier analysis of a shock interacting with three-dimensional isotropic vortical disturbances. This provides the modal structure of the post-shock gas arising from the interaction, along with integral formulas for the amplification of enstrophy, concentration variance, turbulent kinetic energy (TKE), and turbulence intensity across the shock. In addition to confirming known endothermic effects of dissociation and vibrational excitation in decreasing the mean post-shock temperature and velocity, these LIA results indicate that the enstrophy, anisotropy, intensity, and TKE of the fluctuations are much more amplified through the shock than in the thermochemically frozen case. In addition, the turbulent Reynolds number is amplified across the shock at hypersonic Mach numbers in the presence of dissociation and vibrational excitation, as opposed to the attenuation observed in the thermochemically frozen case. These results suggest that turbulence may persist and get augmented across hypersonic shock waves despite the high post-shock temperatures.C.H. was funded by a 2019 Leonardo Grant for Researchers and Cultural Creators awarded by the BBVA Foundation, and by the MICINN Grant No. PID2019–108592RB-C41

Weak-shock interactions with transonic laminar mixing layers of fuels for high-speed propulsion

This paper extends to transonic mixing layers an analysis of Lighthill ("Reflection at a Laminar Boundary Layer of a Weak Steady Disturbance to a Supersonic Stream, Neglecting Viscosity and Heat Conduction," Quarterly Journal of Mechanics and Applied Mathematics, Vol. 54, No. 3, 1950, pp. 303-325.) on the interaction between weak shocks and laminar boundary layers. As in that work, the analysis is carried out under linear-inviscid assumptions for the perturbation field, with streamwise changes of the base flow neglected, as is appropriate given the slenderness of the mixing-layer flow. The steady-disturbance profile is determined by taking a Fourier transform along the longitudinal coordinate. Closed-form analytical functions for the pressure field are derived in the small- and large-wave-number limits, and vorticity disturbances are obtained as functions of the pressure perturbations. The analysis is particularized to ethylene&-air and hydrogen&-air mixing layers, for which the dynamics are of current interest for hypersonic propulsion. The results provide, in particular, the effective distance of upstream influence of the pressure perturbation in the subsonic stream. The resulting value, which scales with the thickness of the subsonic layer, is much smaller than the upstream influence distances encountered in boundary layers. This study may serve as a basis to understand shock-induced autoignition and flameholding phenomena in simplified versions of non-premixed supersonic-combustion problems.This work was supported by the U.S. Air Force Office of Scientific Research grants FA9550-12-1-0138 and FA9550-14-1-0219. We aregrateful to Amable Liñán for useful conversations at the early stages of this project

Theoretical studies in spiral edge-flame propagation and particle hydrodynamics

Applied mathematics techniques are used in this investigation to gain insight into three different physical processes of current interest in combustion and fluid dynamics. The first problem addresses the propagation of spiral edge flames found in von Karman swirling flows induced in rotating porous-disk burners. In this configuration, a porous disk is spun at a constant angular velocity in an otherwise quiescent oxidizing atmosphere. Gaseous methane is injected through the disk pores and burns in a flat diffusion flame adjacent to the disk. Among other flame patterns experimentally found, a stable, rotating spiral flame is observed for sufficiently large rotation velocities and small fuel flow rates as a result of partial extinction of the underlying diffusion flame. The tip of the spiral can undergo a steady rotation for sufficiently large rotational velocities or small fuel flow rates, whereas a meandering tip in an epicycloidal trajectory is observed for smaller rotational velocities and larger fuel flow rates. A formulation of this problem is presented in the equidiffusional and thermodiffusive limits within the framework of one-step chemistry with large activation energies. Conditions for extinction of the underlying uniform diffusion flame are obtained by using activation energy asymptotics. Edge-flame propagation regimes are obtained by scaling analyses of the conservation equations and exemplified by numerical simulations of nearly straight two-dimensional edge flames near a cold porous wall in a von Karman boundary layer, for which lateral heat losses to the disk induce extinction of the trailing diffusion flame but are relatively unimportant in the front region, consistent with the existence of the cooling tail found in the experiments. The propagation dynamics of a steadily rotating spiral edge is studied in the large-core limit, for which the characteristic Markstein length is much smaller than the distance from the center at which the spiral tip is anchored. An asymptotic description of the edge tangential structure is obtained, spiral edge shapes are calculated, and an expression is found that relates the spiral rotational velocity with the rest of the parameters. A quasistatic stability analysis of the edge shows that the edge curvature at extinction in the tip region is responsible for the stable tip anchoring at the core radius. Finally, experimental results are analyzed, and theoretical predictions are tested. The second problem analyzes, in the limit of small Reynolds and ionic Peclet numbers and small clearances, the canonical problem of the forces exerted on a small solid spherical particle undergoing slow translation and rotation in an incompressible fluid moving parallel to an elastic substrate, subject to electric double-layer and van der Waals intermolecular forces, as a representative example of particle gliding and of the idealized swimming dynamics of more complex bodies near soft and sticky surfaces in a physiological solvent. The competition of the hydrodynamic, intermolecular and surface-deformation effects, induces a lift force, and drag-force and drift-force perturbations, which do not scale linearly with the velocities, and produces a non-additivity of the intermolecular effects by reducing the intensity of the repulsive forces and by increasing the intensity of the attractive forces. Reversible and irreversible elastohydrodynamic adhesion regimes are found, and elastohydrodynamic corrections are derived for the critical coagulation concentration of electrolyte predicted by the the Derjaguin-Landau Verwey- Overbeek (DLVO) standard theory of colloid stabilization. The third problem addresses the dynamics of pollen shedding from wind-pollinated plants, and establishes a fluid-dynamical framework for future refinements. A simple scaling analysis, supported by experimental measurements on typical wind-pollinated plant species, is used to estimate the suitability of previous resolutions of this process based on wind-gust aerodynamic models of fungal- spore liberation. According to this scaling analysis, unsteady boundary-layer forces produced by wind gusts are found to be mostly ineffective since the Stokes-Reynolds number is a small parameter for typical anemophilous species and wind streams. A hypothetical model of a stochastic aeroelastic mechanism, initiated by the atmospheric turbulence typical of the micrometeorological conditions in the vicinity of the plant, is proposed to contribute to wind pollinatio

Theory of the propagation dynamics of spiral edges of diffusion flames in von Kármán swirling flows

This analysis addresses the propagation of spiral edge flames found in von Kármán swirling flows induced in rotating porous-disk burners. In this configuration, a porous disk is spun at a constant angular velocity in an otherwise quiescent oxidizing atmosphere. Gaseous methane is injected through the disk pores and burns in a flat diffusion flame adjacent to the disk. Among other flame patterns experimentally found, a stable, rotating spiral flame is observed for sufficiently large rotation velocities and small fuel flow rates as a result of partial extinction of the underlying diffusion flame. The tip of the spiral can undergo a steady rotation for sufficiently large rotational velocities or small fuel flow rates, whereas a meandering tip in an epicycloidal trajectory is observed for smaller rotational velocities and larger fuel flow rates. A formulation of this problem is presented in the equidiffusional and thermodiffusive limits within the framework of one-step chemistry with large activation energies. Edge-flame propagation regimes are obtained by scaling analyses of the conservation equations and exemplified by numerical simulations of straight two-dimensional edge flames near a cold porous wall, for which lateral heat losses to the disk and large strains induce extinction of the trailing diffusion flame but are relatively unimportant in the front region, consistent with the existence of the cooling tail found in the experiments. The propagation dynamics of a steadily rotating spiral edge is studied in the large-core limit, for which the characteristic Markstein length is much smaller than the distance from the center at which the spiral tip is anchored. An asymptotic description of the edge tangential structure is obtained, spiral edge shapes are calculated, and an expression is found that relates the spiral rotational velocity to the rest of the parameters. A quasiestatic stability analysis of the edge shows that the edge curvature at extinction in the tip region is responsible for the stable tip anchoring at the core radius. Finally, experimental results are analyzed, and theoretical predictions are tested